Magnetic field imaging with atomic Rb vapor 
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We demonstrate the possibility of dynamic imaging of magnetic fields using electromagnetically 
induced transparency in an atomic gas. As an experimental demonstration we employ an atomic 
Rb gas confined in a glass cell to image the transverse magnetic field created by a long straight wire. 
In this arrangement, which clearly reveals the essential effect, the field of view is about 2x2 mm 2 
and the field detection uncertainty is 0.14 mG per 10 /im x 10 fim image pixel. 
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Efforts to develop measurement techniques for preci- 
sion magnetometry have a long and successful history 
in atomic physics [J, 0, H, 0, [j| . Some recent advances 
in development of compact atomic magnetometers take 
advantage of the unique properties of the electromagneti- 
cally induced transparency (EIT) effect and the existence 
of coherently produced optical dark states, which can be 
extremely sensitive to external fields @. Here, we pro- 
pose the idea of a three axis magnetometer for dynamic 
imaging of magnetic field gradients based on EIT using 
state of the art imaging technologies. We also demon- 
strate the essential idea with a one axis (two dimensional) 
prototype. As described later, a fully implemented three 
dimensional magnetic gradiometer promises both high 
precision and useful visualization tools for both static 
and dynamic magnetically complex environments. The 
potential applications of this technique range from es- 
sentially fundamental (e.g. precision searches for mag- 
netic monopoles 0) to quite applied, such as, for ex- 
ample, precise monitoring of the magnetic environment 
in medical applications or quantum information experi- 
ments Hi]. 

The proposed approach relies on the spectrally ultra- 
narrow transmission resonances observed under EIT con- 
ditions in an atomic gas. These resonances are associ- 
ated with the coherent non-interacting (dark) states in a 
three-level A system formed under the combined action 
of two electromagnetic fields in two-photon Raman reso- 
nance with two long-lived metastable states of an atom. 
The width of the EIT resonances is determined by the in- 
tensities of the laser field and the lifetime of the ground- 
state coherence, and can be as small as a few Hz jlSEI- 
Such spectral sensitivity makes this effect very promis- 
ing for precision metrology. In particular, the spectral 
location of the transparency resonance depends on the 
external magnetic field when the A-system is based on 
the magnetic field-sensitive Zeeman sublevels. Moreover, 
if the magnetic field is spatially varying transversely in a 
plane orthogonal to the probe beam's k vector, the res- 
onances will occur at different spatial locations for dif- 
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FIG. 1: A schematic diagram of the experimental configura- 
tion. 



ferent two-photon resonance detunings. Thus, the spa- 
tial distribution of transmitted intensity will display a 
transverse image of the spatial regions where the two 
photon resonance condition is satisfied. Recording a se- 
ries of such images for various two-photon frequencies 
can create a spatial map of the transverse magnetic field. 
Combining such measurements made in three orthogonal 
directions (formed from, for example, three independent 
arrangements) allows six elements of the VB tensor to be 
measured. According to the source free Maxwell's Equa- 
tions, there are four constraints on the nine components 
of VB, implying a desirable unit redundancy in each set 
of six measurements. 

In this Letter we illustrate the concept by imaging the 
magnetic field produced by a long straight wire carrying 
a steady current / running along the length of a Rb va- 
por cell. Since in our experiment we are able to detect 
the light transmission through the cell only in one direc- 
tion, we took extreme care to make the wire as coaxial 
as possible with the laser beam to remove any variation 
of the magnetic field along the optical axis. The aim 
is then to measure the resulting transverse variation of 
the magnetic field across the laser beam by imaging the 
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magnetic field across the laser beam is mainly along the 
y-axis and changes linearly with small displacement 5x: 



m=-2 m=-1 m=0 m=1 m=2 



FIG. 2: Energy level diagram illustrating the Zeeman transi- 
tions essential to the EIT A schemes formed by the carrier Eo 
and the first modulation sideband E+i. 



transverse variations of the transmitted light intensity as 
a function of the two-photon detuning. 

A schematic of the experimental setup is shown in 
Fig. [T] We use a vertical-cavity surface emitting laser 
(VCSEL) directly phase-modulated at 6.834 GHz (a 87 Rb 
hyperfine frequency), and use the fundamental (carrier) 
laser frequency and one of the side bands to produce 
the two optical fields required for EIT observation. The 
laser beam with total power 200 pW and a slightly el- 
liptical Gaussian profile [waist sizes 1.8 mm and 1.4 mm 
full width half maximum (FWHM)], was directed into 
the cylindrical Pyrex cell (with length 75 mm and diam- 
eter 22 mm) containing isotopically enriched 87 Rb vapor 
and 15 Torr of Ne buffer gas. The cell was mounted in- 
side a three-layer magnetic shielding to reduce stray lab- 
oratory magnetic fields with the suppression factor of at 
least 1000, and its temperature was actively stabilized at 
322 K by regulating the current through a bifiler heater 
wire wrapped around the innermost shield. 

The magnetic field producing wire was mounted par- 
allel to the laser beam at a distance po = 20.1 ± 0.2 mm, 
and connected to a low-noise current supply operating at 
I = 438 ± 1 mA. Then the magnetic field outside the wire 
produced by this current is orthogonal to the light propa- 
gation direction, and can be estimated as (neglecting the 
effects of the magnetic shielding): 
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where p is the distance from the center of the wire, and 
po is permeability of free space. If we let the laser beam's 
wave vector and the direction of current flow define the 
z-direction, the curves of constant magnetic field B are 
then circles in the x-y plane and centered on the wire. 
For a given two-photon detuning the EIT resonance con- 
ditions are obeyed only for a given value of the magnetic 
field, and thus for a large laser beam we in general ex- 
pect to observe a corresponding "bright" circular arc in 
the transmission spatial profile, with the radius of the arc 
depending on the set two-photon detuning. In our exper- 
imental arrangements, however, the laser beam diameter 
is small compare to the distance to the wire po, so the 
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For the measurements, the laser frequency was tuned 
such that the stronger optical field (unmodulated carrier) 
was resonant with the 5 Si/ 2 F — 2 — > 5 P%/2 F 1 — 1 
transition of 87 Rb, while the frequency of the +1 mod- 
ulation sideband matched the 5 S\/2 F = 1 — > 5 P\ji 
F' = 1 transition. The laser output was linearly po- 
larized. In this configuration there are two possible A- 
systems for F — 1, 2; nip = ±1 — > F' = 1; mp' = ±1 that 
produce EIT resonances at two-photon detunings Av — 
A HFS ±2g m B (shown in Fig.H]), where g m = 0.7 MHz/G 
is the gyromagnetic ratio for 87 Rb [HI]. We also adjusted 
the power of the VCSEL modulation signal to cancel the 
light shift of the EIT resonances, and to avoid the effects 
of non-uniform spatial distribution of the laser beam on 
the measured resonance position [l4l ]. For our experi- 
mental conditions such cancelation occurred for a side- 
band/carrier intensity ratio of 0.75. Also, we analyzed 
only the area of the image where the intensity of the 
beam exceeds half of its maximum value. 

To map the spatial variation of the magnetic field, we 
stepped the VCSEL modulations frequency in a 20 kHz 
range around one of the EIT resonances in 200 Hz in- 
crements, taking images of the transmitted laser beam 
profile using a digital camera (Unibrain Fire-i 511b). 
The dominant source of noise in our experiment was due 
to the phase-to-amplitude noise conversion of the VC- 
SEL large phase noise (the lincwidth of the laser was w 
100 MHz) by the absorption in the cell. To reduce the 
detection noise of the image, we recorded and averaged 
200 sequential images grabbed from the camera, which 
operated at a 30 frame per second rate for any particular 
two photon detuning. Even then the detection noise is 
larger than the digital resolution of the 12 bit analog- 
to-digital converter of the camera. For each pixel of the 
image we plotted the intensity as a function of the two- 
photon detuning, as shown in Fig. [3](a) ■ Then we located 
the position of the EIT maximum (and correspondingly 
the value of the local magnetic field). After analyzing 
the transmission for every pixel we obtained a spatial 
map of magnetic field, shown in Fig. [2Kb). The average 
variance between two sequential runs allowed us to put 
an upper limit on the experimental uncertainty of the 
magnetic field measurements AB = 0.14 mG for every 
10 pin pixel. 

We compare the measured magnetic field with numer- 
ical calculation of the magnetic field generated by a long 
straight wire inside a cylindrical magnetic shield. Since 
the inner shield modifies the magnetic field from the sim- 
ple Biot-Savart law given by Eq. [TJ we used a finite ele- 
ment (FEM) analysis program to calculate the expected 
magnetic field distribution inside the laser beam [la ], as 
shown in Fig. [3Jc). The rather large uncertainty of our 
FEM model is governed by possible variations in each 
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FIG. 3: (a) Transmission through the cell detected at three different camera pixels [marked in (b)] as functions of the two- 
photon detuning. Solid lines show low-pass filtered data used to determine the maximum transmission position, (b) Map of the 
magnetic field strength. Dashed lines mark the 75% and 50% levels of the maximum laser beam intensity, (c) Experimentally 
measured magnetic field slice of the map (b) along y = as a function of x position, and the theoretically calculated magnetic 
fields using the Biot-Savart law [Eq.Q] and a finite element analysis. The error bars in the FEM curve represent the variation 
in the calculated magnetic field due to ±1 mm uncertainty in the magnetic shielding position with respect to the wire and the 
laser beam. 



shield's diameter, thickness, and imperfect cylindrical 
shapes, as well as the uncertainty in the location of the 
wire inside of the shielding. For instance, Fig. [3]Jc) illus- 
trates that displacement of the shielding position as small 
as 1 mm is sufficient to match the FEM simulations with 
the experimental measurements. 

The demonstrated sensitivity of our apparatus (~ 
lAfiG/fim) is limited by the transmitted intensity fluc- 
tuations of the broadband VCSEL output. We expect 
much better sensitivity for narrow-band electromagnetic 
fields, which should dramatically reduce the amount of 
amplitude noise in the output. A good candidate for such 
a source, for example, could be a narrow band external 
cavity diode laser followed by an external phase modula- 
tor. 

We can estimate the precision of the magnetic field 
measurements AB due to the uncertainty in two- 
photon detuning frequency corresponding to the maxi- 
mum transmission for a given signal-to noise ratio S/N 
0: 

In this expression we assume that the EIT resonance has 
Lorentzian shape with a known FWHM 7, and n is the 
number of sampling points per resonance width. The 
optimal frequency span in this case should exceed, or be 
comparable with, the value 7, with n > 5 — 10 for each 
recorded trace to avoid numerical errors in fitting related 
to discretization artifacts. 

We note that the ultimate relative sensitivity of our 
method should be similar to that of atomic clocks, based 
on the same operational principle [12j |. In particular, 
miniature atomic clocks demonstrate a fractional stabil- 



ity on the order of 10~ n — 1CP 12 / \/Hz. However, in prac- 
tice the performance of the magnetic field imager may be 
limited by other factors, such as digital noise or the lim- 
ited optical sensitivity per pixel of a charge coupled cam- 
era. These, and other factors including the essential role 
of longitudinally varying magnetic fields, are currently 
under investigation. 

In conclusion, we propose a new method of dynamic 
imaging of a magnetic field by detecting spatial and 
temporal variations in transmission of light through an 
atomic sample under conditions of electromagnetically 
induced transparency. In general, this method can be 
applied to measurements of three-dimensional magnetic 
field maps. To demonstrate the basic operational princi- 
ple, we reconstructed the spatial distribution of the trans- 
verse magnetic field created by a current-carrying long 
straight wire, running along an atomic Rb vapor cell, 
with the precision of 14^G//j,m inside the 1.8 x 1.4 mm 2 
laser beam cross-section. Our prototype apparatus can 
be readily scaled up for either a larger or smaller de- 
tection area by adding a zooming telescope in front of 
the camera (although the total power of the laser would 
have to be reoptimized). Also, since our prototype ex- 
periment was performed using a VCSEL, we expect that 
this method can be used with recently developed chip- 
scale systems for miniature atomic clocks and magne- 
tometers [l7| . 
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